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Motivation for Scaled Experiments 

Brief Scaling History 

- Steady-State Combustion 

- Combustion Stability 

- Life Prediction 

Scaling Approaches Presently Used at Purdue 



Background 

Stated goals for long-life LRE’s have been between 100 and 500 cycles 

- Inherent technical difficulty of accurately defining the transient and steady 
state thermochemical environments and structural response (strain) 

- Limited statistical basis on failure mechanisms and effects of design and 
operational variability 

- Very high test costs and budget-driven need to protect test hardware 
(aversion to test-to-failure) 

Ambitious goals will require development of new databases 

- Advanced materials, e.g., tailored composites with virtually unlimited 
property variations 

- Innovative functional designs to exploit full capabilities of advanced 
materials 

- Different cycles/operations 

Subscale testing is one way to address technical and budget challenges 

- Prototype subscale combustors exposed to controlled simulated conditions 

- Complementary to conventional laboratory specimen database development 

- Instrumented with sensors to measure thermostructural response 

- Coupled with analysis 



SSME Film Cooling Analysis 


Configuration • 

- Propellant = LOX + LH2 with O/F = 

6.02 

- M_dot_LOX = 64,000 liter/min 

- M_dot_LH2 = 178,000 liter/min 

- M_dot_coolant for regen cooling = 

29.06 lb/sec 

Chamber condition 

- Pc = 3300 psi 

- Tc = 3500 K (5840 F) 

- D_throat = 10.88” 

- E = 77 

Cooling channel 

- Wall thickness = 0.03” 

- Width = 0.04 “ 

- Height = 0.12“ 

- Pressure_throat = 3851 psi 

Thermal condition at throat 

- Heat flux = 80 Btu/in A 2-s 

- hg = 58000 W/m A 2-K 

- Twg = 1100 F 

Wall adiabatic temperature 

- Taw = Tr - r|(Tr-Tco) 

Where Tr = recovery temperature 
r) = film cooling efficiency 
Tco = initial coolant temperature 


Current near wall O/F ratio 

- q_dot = hg(Taw-Twg) 

Where q_dot = 80 Btu/in A 2-s 
hg = 58000 W/m A 2-K 
Twg = 1100F 
Taw = 3125 K 
r| = 0.5 
Tco = 2750 K 

4 0/F_nw = 3.54 from Flame temperature 
O/F ratio chart 



SSME Film Cooling Analysis 

SSMEO/Fvs Life 


Current film cooling 
condition 

- 0/F_nw = 3.54 

Parametric study with 
fixed film flow rate (5 %) 

*Porowski et al. method 
(AIAA Journal Vol. 2 No. 
2, 1985) 

- 0/F_nw change = 3.54 — » 

1.0 

- Life change = 61 — > 107 
(75.4% increase) 

- Isp change = 465 — ► 457 
(1.83 % decrease) 



Isp vs 0/F Variation 
(coolant m_dot = 5.0 %) 



Scaling Objectives and Approaches 

Combustor characterization is goal 

- Validation data for design analysis models 

- Assess innovative functional design, materials, operation 

- Investigations into specific physics 

Single element, multi-element, 40K, 250K 

Cold flow and hot fire 

Performance, heat transfer, life, stability 

Experimental objective needs to define scaling 
approach and measurement 

- Well-instrumented combustors linked to analysis 

- Thrust level and number of elements 

- Element scaling and configuration 



difficulty / cost 


Hierarchy of injector experiments 



jector design 


Single element 
atmospheric 
cold flow 


Single element 
high pressure 
cold flow 




Subscale 
Multi-element 
hot fire 






AFRL cold flow facility 


degree of simulation 







Brief History of Scaling in the US - 
Steady State Combustion 


• JPL studies of mixing efficiencies of impinging jets 

• Bell Aerospace/ AFRL holographic and shadowgraphic 
studies of combusting flows 

• Rocketdyne development of LISP methodology for SDER 

• Aerometrics development of PDPA 

• Rocketdyne studies of flameholding behind LOX post 

• PSU measurements of chemical species in HO combustors 

• AFRL studies of supercritical jets 



Single Element Test Chamber 



Stability Scaling 

Simulation of chamber dynamics in subscale 
configuration is very difficult 

- Acoustic frequencies scale as ~ 1/d 

- Pressure v velocity sensitivity 

Scaling approaches 

- Wedges, T-bumers, 2-d chambers 

- IT = 3T scaling 

Single element rarely used in US, but is more 
typical in Russia 




Experimental Approach of Bazarov 


This facility screened 
Injector elements for 
Liq/liq and gas/liq 
Injectors for over 20 
Years (1965-85) 

Typical Pc = 750 psi, 
Total flowrate of 5 lb/s 

‘self-oscillation’ and 
response to pulsations 
measured 



Fig.X Pneumatic and hydraulic scheme of lire stand 

I -piston, 2-measuring vessel, 3-nozzle collector, 4-combustion chamber, 5-pulsator, 6-fuel tank, 7-time 
delay valve, X-blow through valve, 9-main bi-propellant valve. 10,11 -pressurising gas reductors, 1 2- 
oxidizer tank, 13-exhaust tubes, 14-water valve, 15-oxidizer return tank, 16-pressurised water tank, 17- 
ejector, lfi-;nr compressor 
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Use full-scale injector 
elements 

Experiment designed to 

simulate controlling processcooikg Water | 

- mixing 

Match equivalence ratio and 
volumetric flowrates using 
diluted gaseous propellants 

Combustor acoustics 
matched by using 
appropriately sized low- 
pressure chamber 

Stability boundaries 
determined by varying 
flowrates 

Relative boundaries indicate 
stability ranking 


Experimental Approach of 
NIICHIMMASH 

Cooling Waicr 
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L - Combustion Qi amber 

2 * Steel Fiat Plate 

3 - Injector Element with Oxidizer 

and Fuel Manifolds 

4 - Coil 

5 - Muffle Furnace 
- Throtliing Cock 

7 - Metering Orifice 
S - Pulsation Transducer 

9 - Amplifier 

10 - Oscilloscope 

II'- Spectrum Analyzer 

12 * Voltmeter 

13 - Frequency Meter 


Figure 6. Schematic of Single Element Model Set- up and Instrumentation 
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Propellant Distribution Effects 
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Table 7 - 2 . 5 a. — G as Rocket Test History With 
VART 6 U 3 Injection Profiles 


Instabilities initiated apiaiitiuuioLisly unci lircarlv; mean 
clumber pressure, 150 psia; combustion clumber diameter, 
7 in.; combustion chamber length, r t In.] 
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Figure 7 .^ 5 m- — Injection rndrn! profile cdiiipdriSOIfE 





Single Element ‘Instability’ 



Combustor oscillations at 
driven atomization frequency 


Impinging jets driven by 
piezoelectric actuator 



Time (*) 




Subscale Test Activities at Purdue - 
Maurice Zucrow Laboratory 





Advanced Propellants and 
Combustion Lab 



Two cells w/ 1 Klbf thrust stands 
Propellant supply of 1800 psia 
2-4 gallon oxidizer tanks 
1 & 4 gallon fuel tanks 

National Instruments hardware & 

Lab View software 

- 32 channels pressure 

- 32 channels temperature 


•All valves computer 
controlled 

•Rapid test article 
installation 

•Design/Build/Test course 



High Pressure Lab 

Renovation funded thru Indiana 21 st Century R & T Fund - 
Propulsion and Power Center of Excellence 
Facility activated in May ‘03 













6,000 psi Nitrogen System 



• Pressurization, Actuation and Purge Gas 

• 2,400 gallon Liquid Nitrogen Tank w/ 6,000 psi Pump 

• 253 ft 3 6,000 psi Nitrogen Tube Trailer 

• Computer Controlled Pressurization Systems 


Propellant/Coolant Tanks 



• 22 gal 5,000 psi LOx 

• 16 gal 5,000 psi Fuel 

• 220 gal 5,000 psi H 2 0 

• 400 gal 800 psi H 2 0 2 

• Hydraulic Control Valves 





10,000 lbf Thrust Test Cell 


• LabView 6.1 -based DACS 

• 10,000 lbf thrust measurement 

• 64 channels pressure 

• 96 channels thermocouples 

• 18 channels analog control 

• 32 channels on/off control 




Control System Operation 


• Data System Located Adjacent to Test Cell 

• Operation Remoted to Control Room (KVM Extender) for Testing 

• Video Recorded Directly to DVD 




Test Cells 



18” Thick Reinforced Concrete Test Cell Walls 
High Flow Capacity Test Cell Exhaust Fans 
Heated High Pressure Air Plumbed to Both Cells 
Walled Containment Area 


Injector Characterization Scaling Approach 

•Study Objectives 

-Steady state and dynamic characterization of ORSC MC 
injector elements 

•Approach 

-Investigate full-scale elements at realistic operating 
conditions 

-No film cooling (if possible) 

-Evaluate different injector design configurations 
-Couple with analysis 

•Measurements 

-Energy release profile from axial pressure gradient 
-Injector face and chamber wall thermal environments 
-Plume signature with IR tomography 
-Manifold, injector and chamber p’ 



ORSC Main Combustor Components 
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271 elements, 1722 lb f each, d = 0.5 in 



Principle Design Features 






Gasflow inlet lip to affect 
acoustic admittance 


Swirl injectors isolated 
from gas path 


Ox tube length set to 
tune injector acoustics 


Beveled recess to 

provide injector face 
cooling 


Two rows of inlets 

provides wave 
cancellations 


Atomization due to 
Kelvin-Helmholtz instability 

Two-phase region enhances 
impedance 



Liquid submergence 

to enhance mixing, 

control face heating & 
promote stability 
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Single Element Sizing Exercise 



Approach 

• use full scale F/element (1722 lb fvac ) 
mox = 3.6 lb/s, mf = 1.2 lb/s 


test at ‘full’ Pc (2250 psia) 

At = 0.39 in 2 , dt = 0.70 in 
match injection pressure drops (10%) 
czj> dinj, ox = 0.43 in, dinj = 0.57 in 


Possible scaling methods: 
Contraction ratio (1.61) i 


f dc = 0.89 in 


Element to chamber area ratio (0.30) cz) dc = 1.04 in 
Element-element spacing (0.60d) i=) dc = 0.9 1 in 
Element- wall spacing (0.60d ?) m) dc = 0.91 in 
Element area (0.65 in 2 ) i=J> dc = 0.91 in 
Chamber length based on L* ~ 30 in (??) 
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Baseline Injector Design 


MATING FLANGE / 

Sensor Ring RP-I Inlet 

\ ^ 



InLet / X 

P' Measurement Injector Face 



High-Pressure 

Chamber 


Igniter Section 

Calorimeter 

Sections 

Flanges 




Jacket 



Gun-Drilled 
Chamber Sections 


Calorimeter 

Sections 


Nozzle 

Liner 


Gun-Drilled 
Chamber Sections 



Dump-Cooled 
Throat Section 


Life Prediction 

Rocket combustor liner such as SSME 
operated at high temperature (6000F) 
and pressure (3000 psi) ranges as well 
as extreme heat flux (80 Btu/in 2 -s) 
requires active cooling devices to 
prevent material failure. 

Combustor liner experiences high 
thermal structural stress (-100 MPa) 
during mission profile (SSME 8 min) 

Experiments by Quentmeyer and 
Jankovsky showed bulging and 
thinning of liner due to cyclic loading 

Kasper and Porowski developed 
analytical life prediction methods 
using simple fatigue and creep model 

Robinson, Arnold and Freed 
developed visco-plastic model for 
fatigue-creep interaction phenomena 
which is believed to be a main failure 
mechanism 


Background 



Typical failure mode of combustor 
liner at throat so called “dog house 
effect” per Quentmeyer 






Full Scale - Subscale Life Comparison 

- Pc = 3300 psi, Tc = 3500 K 





Full scale engine 

Strain_max = 2.4 
Life = 120 



Strain_max = 3.94 
Life = 48 










Approach 

Develop DBT course with life prediction as part of AAE curriculum 
Develop design requirements 

- Controlled hot-gas environments - use ‘pre-combustor’ 

- Creep-fatigue interaction failure of cooled liner 

- Failure within reasonable number of cycles 

Life prediction analysis using conventional methods 

- Chemical equilibrium in pre-combustor 

- One-dimensional heat transfer analysis for initial design 

• critical heat flux and cooling requirements, duty cycle 

- FEM for stress and plastic strain 

- Strain-life curves for cycle life 

- More advanced life modeling by graduate student following project 
Cyclic testing of test article 

- Ten cycles per test 

- Validation of cooling analysis 

- Regular inspection 
Test-to-failure 



Combustor Design Parameters 


• Top level requirements 

- Less than 200 life cycle 

- Test should produce verifiable 
results 

- Liner has no melting prior to the 
LCF failure 

- All parts had to be manufactured 
in ASL at Purdue 

• Under these requirements, the 
coolant pressure, flow rate and 
cooling channel aspect ratio (0.5) 
were determined. 



Parameter 

Value 

Propellant 

90% H 2 0 2 + JP-8 

Propellant mixture ratio (O/F) 

4.0 

Propellant flow rate 

1.25 lb/s 

Chamber pressure (P c ) 

200 psia 

Chamber temperature (T c ) 

3440 °F 

Characteristic velocity (C*) 

4961 ft/s 

Throat area (A t ) 

0.915 in 2 

Characteristic length (L*) 

70 

Test liner diameter 

2.0 in 

Test liner length 

5.0 in 

No. of cooling channel 

30 

P 

coolant 

110 psi 

M - dOt C00l», 

0.8 lb/s 


Table 1 : Combustor design parameters 



Twg & Twl [K] q [Btu/in A 2-s] 


Thermal Structural Prediction 


Thermal analysis 

• Burn out heat flux — 6.54 Btu/in 2 

• Max wall temp — 670 K 





Axis [in] 



Total strain 
predicted by 
ANSYS around 
rectangular cooling 
channel. 

-Total strain — 2.0 % 
-Life expected — 115 
cycles 


Strain-life curve for 
OFHC at 810 K 
from NASA CR- 
134806, 1975 



Test Article 



Catalyst bed for decomposing H 2 0 2 

Heat sink dump combustor for hot gas generation 

Chamber liner — water cooling 

Center body — water cooling with TBC (0.01” 
thick) 


To Ox Main Valve 
500 psi 90% H202 
Vi ” AN Fitting 


P_Oxcatout 
14” AN Fitting 
200 psi 

Engine Mount Bolts to 
(4) Unitstrut L-brackets 
On Test Stand w/(4) Vz bolts 


T_Precombustor 

Welded 

260 degrees F 

P_Chamber1 
14" AN Fitting 
200 psi 


P_CBin 
14" AN Fitting 
100 psi 

T_Jackin 

1/16” Swagelok Fitting 
71 degrees F 


P_Jackout 
14” AN Fitting 
80-100 psi 



T_Catout 

1/16” Swagelok Fitting 
1200 degrees F 


To Fuel Main Valve 
250 psi RP-1 
1/2” AN Fitting 


P_Chamber2 
14” AN Fitting 
200 psi 


To Water Main Valve 
200 psi H20 
1/2” AN Fitting 


P_Jackin 
14” AN Fitting 
100 psi 

T_Jackout1 

1/16” Swagelok Fitting 

150 degrees F 

T_Jackout2 

1/16” Swagelok Fitting 

150 degrees F 




Testing 

• Tests were conducted in the APCL at Purdue University 

• Propellant flow timing sequence was automatically controlled by 
pneumatically actuated valve with LAB VIEW system 




Test article assembly on test stand 


Cyclic test 



Test Results 


Chamber pressure, C* efficiency, propellant 
mass flow rate, coolant temperature and 
pressure were measured and calculated 

Data reduction was performed using in-house 
code written by students using MATLAB 

Validation procedure 

- Measure coolant AT, wall thinning rate 

• 2.15E-5 in/cycle (0.032”^0.029”) 

- Verify ID thermal model 

- Compute updated thermo- structural 
environment 

- Make life prediction 


Casing Jacket T errsps. Life Cycle Anatysis: Cycle T esling 1 2. S May 2003 



Coolant temperature 



Predicted and measured coolant temperature 
AT = 4.0K at throat 



Discoloration and deformation at 


90 cycles (1. 5” X 0.6”) 




Updated Structural Analysis 




Simulation of temperature, strain and deformation (bulging, thinning) using 
ABAQUS explicit module 

Maximum strain : 1.2 % at middle of ligament 
Only bulging of ligament was simulated 


Deformation after 80 cycle 


Deformation after 100 cycle 


Deformation after 60 cycle 







Summary and Conclusions 


Small-scale rocket combustor was 
designed and tested to verify life 
prediction models for low cycle 
fatigue and fatigue-creep interaction. 

Several life prediction methods were 
applied to predict combustor life and 
were compared with test results. 

Correlation data used to improve 
predictions. 

Improvements would include fixing 
the liner lands to the structural jacket, 
and testing at more severe conditions. 


Prediction 

method 

Estimated life 
cycle 

Determined 
life cycle by 
experiment 

Effective 

stress-strain 

115 


ANSYS 

115 

270 

Porowski 

51 

Dai and Ray 
with Freed 
model 

260 


ABAQUS 

320 


Comparison of life prediction with test 


Summary and Conclusions 


100’s of cycle goal is very challenging and 
verification would be very expensive 

- Question of economic feasibility 

Improved life prediction methodology for 
expanding range of design and operational 
scenarios is needed 

- Probabilistic life prediction design analysis 

- Testing methodologies with in situ thermo structural 
response measurements 

- Environments definition 

- Improved material database and understanding of 
damage mechanisms 
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